Preferred Measurements: Optimality and Stability in Quantum Parameter Estimation 



Gabriel A. Durkin' Q 

'Quantum Laboratory , NASA Ames Research Center, Moffett Field, California 94035, USA 

(Dated: September 11, 2009) 

We explore precision in a measurement process incorporating pure probe states, unitary dynamics and com- 
plete measurements via a simple formalism. The concept of 'information complement' is introduced. It under- 
mines measurement precision and its minimization reveals the system properties at an optimal point. Maximally 
precise measurements can exhibit independence from the true value of the estimated parameter, but demanding 
this severely restricts the type of viable probe and dynamics, including the requirement that the Hamiltonian 
be block-diagonal in a basis of preferred measurements. The curvature of the information complement near a 
globally optimal point provides a new quantification of measurement stability. 
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Scientists strive for an understanding of Nature by a physi- 
cal interaction introducing correlations between observer and 
observed, and the process is called measurement. Fundamen- 
tal limitations on the precision of any measurement exist due 
to the geometric distinguishability of quantum states (T), yet 
for a quantum observable there is still a classical probability 
distribution over measurement outcomes. It may depend on 
some real-valued system parameter 9 such as an interaction 
time or interferometric phase that has no associated Hermitian 
observable and is not measurable directly. Yet, the estimation 
of 9 may be the true goal of the measurement. Inferring 9 
within a confidence interval (precision) A8 from frequencies 
of measurement outcomes is a standard challenge in classical 
information theory with an established methodology Q. In 
a quantum context, the conventional approach stretches back 
decades [3|, employing techniques from Riemannian geome- 
try to find precision limits via 'quantum Fisher information' 
(QFI) [4, 5j|6][7] , a function of the input state and dynamics 
alone. It defines a precision limit for measurements but sheds 
little light on what measurement to use - those proposed as 
optimal are typically functions of the unknown 9 [8|. (It is 
of limited utility that the parameter estimation should require 
prior knowledge of its true value. Adaptive techniques em- 
ploying feedback have been proposed to circumvent this prob- 
lem 0.) 

In this letter we harness a quantum formalism for pure 
states emerging directly and naturally from a single result 
of classical information theory. The method incorporates all 
three instrument components - input, dynamics (Hamiltonian) 
and measurement choice - on an equal footing without re- 
course to QFI or the mathematical apparatus supporting it. 
Important metrological results will be confirmed in a straight- 
forward manner along the way towards developing new con- 
ditions for optimality and stability in quantum measurements. 

A preferred measurement M has three important proper- 
ties: (1) high precision, or in the many-body case, 'supra- 
classical' precision. (By this it is meant that the precision is 
better for collective dynamics than is possible for component 
parts evolving in a separable state J6).) (2) If possible, M is 
independent of the estimated parameters |5][T0), and (3) M 
is highly stable, i.e. the precision exhibits robustness against 



small perturbations in the state or measurement alignment. 
Below, these goals are given explicit mathematical expression; 
optimality criteria imposed on probe and dynamics. 

Consider a maximal test [ 1 1 1 having outcomes labelled k 
and an associated probability distribution P(9) — {pk(8)} 
that depends on a continuous real parameter 9. The parame- 
ter estimation task involves an inference of 9 from {pk(9)}. 
Classical Fisher information 1 12] is defined as 

J (9) = J2Pk(0)[d s In Pk (9)] 2 , (1) 

k 

a measure of the information contained in the distribution 
P(9) about the parameter 8 0. Unlike the quantum coun- 
terpart [1| it defines a unique distance metric on probability 
space [13]. An explicit lower bound for the standard error of 
an unbiased estimate 9 on the true value 9 is given by the re- 
ciprocal of the Fisher information, (S9) 2 > \/J(9), called 
the Cramer-Rao bound lfl"2"l . For optimal precision one must 
therefore maximize J (9). 

Take a complete measurement observable M with out- 
comes {m k } associated with distribution {p k }. Apply M to 
a quantum system that previously evolved from a known ini- 
tial state \ipo) under the dynamics of some time-independent 
Hamiltonian H for time 9. The Schrodinger equation gov- 
erns the dynamics idg\ipe) = H\%jjg), (where li — 1) and 
the time evolution is explicitly \ipg) — exp{— iH9}\ipo) . 
Writing the spectral decomposition of the measurement as 
M = 'Yl lk 'm k \k)(k\ then complex amplitudes (k\ipe) — 
r k exp{i(j) k } give probabilities p k = (k\ip e ) (ipe\k) = r\ , 
where {pk,fk} € [0, 1] and <fi k <E [0,27r) are all real-valued 
functions of 8. Replacing p k with r\ in Eq.([T} gives: 

J(8)=AY,r 2 k , (2) 

k 

which we can now use to find an operator expression for the 
classical Fisher information. Differentiating r k gives 

2r k r k = (k\^ e }(^ e \k) + (A# e )(^|fc) . (3) 

Now, (k\ip e ) = d e {r k e i<t,k ) = e^(r + ir<f> k ) = 
| (k\ipg) |e 4 ^ fc+rfc ) where we define a velocity vector with ra- 
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FIG. 1 : An optimal measurement scheme for a Hamiltonian H span- 
ning any number of dimensions may be restricted to the qubit sub- 
space of its extremal eigenvectors |Af) and | Aj.) . (Other measure- 
ment elements G {\k}} span an orthogonal subspace.) On the Bloch 
sphere: the diagonal basis of H defines the z-axis, and the opti- 
mal measurement projectors \k±) the x-axis. An optimal probe state 
l^opt) lies anywhere in the equatorial plane. 



where we define the 'information complement' : 

no) = E( r ^) 2 = = 



(8) 



a non-negative functional of the probe, measurement and 
Hamiltonian that can only reduce J . Here <p is a classical 
random variable taking values from the set {<j>k} and the c 
subscript denotes the expectation value is classical. We now 
look for a basis {\k)} and associated set {r k , <f>k} that mini- 
mizes JC for a fixed input \ipo)- First, a new description of the 
Hamiltonian expectation value is needed: 



dial and transverse components {fk, r^k], and inclination: 

r fc = tan _1 (r fc fe /f fe ) = axg(k\ip e ) - a,rg(k\ipg) . (4) 

Eq.|3]l yields f k = cos Tk\(k\ipg)\. From the 
Schrodinger equation we also have i(k\ipg) — (k\H\ipg). 
Substituting this expression and squaring gives r 2 = 
cos 2 Tk(ipg\H \k) (k\H \ipg). Summing over all outcomes 'k' 
gives the Fisher information using Eq.(|2|: J{&) = 
^Sfc cos2 T k{^e\H\k) (k\H\ipg). Thus we can define a non- 
linear, positive and hermitian operator Fg, diagonal in the 
measurement basis: 

F e = 4j]cos 2 T k \k)(k\ = Aj2^k\k){k\ , (5) 

k k 

such that the Fisher information is then 



J{6) = {ll>g\HFgH\ll>g) = {^o\H^gH\^ } 



(6) 



where <f> 9 = e lH6 F e e' lH0 = ^Y,k c ^A k '){k'\ is the uni- 
tarily transformed operator. Due to the non-linear nature of 
Fg, a basis transformation \k) i— > \k') gives a different result: 
Fg >—> Fg = i'^2 k c^ )) k'\k'){k'\ , not equivalent to Qg, since 
c-ip.k' 7^ Ctl>,k generally. This explicit definition of Fg will be 
useful in the qubit optimization to come. 

Fixed Probe Optimization: Now let us establish an upper 
bound for J(9) in terms of a fixed input and dynamics, but 
varying the measurement. The completeness of the measure- 
ment basis provides a resolution of the identity J2k\k)(k\ = 
1, and therefore, with cos 2 = 1 — sin 2 r k , Eq.([6]l becomes: 



(H 2 ) -5>in 2 T k (yj g \H\k)(k\H\<p g ) 

k 

(H 2 )-J2(rk<j>k) 2 = (H 2 )-IC(e), 



(7) 



■(<P)c 



(9) 



Therefore, JC - (H) 2 = A 2 ^ > 0. Comparing Eq.0 it fol- 
lows directly that J < 4A 2 H for a fixed input \ipo)- (This 
was derived by a different method in |5|.) The bound is satu- 
rated by a particular qubit input, as we will now show. 

Optimizing for a Qubit: Probes that are eigenstate of H 
give fk i— > and J = from Eq.|2) because any eigen- 
state | A) of H only gains a phase during its evolution. Thus 
r k (6) = |(fc|e^ e |A)| = |(fc|e^|A)| = |(fc|A)| = r*(0), 
and fk = 0. For optimality over all \^>q) and {\k)} the input 
state must thus be a superposition of at least two Hamiltonian 
eigenvectors, |^o) l— * COS7IA1) + e lx sin7|A2). Choose such 
a qubit probe and a measurement basis ^2)} spanning 

the same (D 2 as |Ai,2) : \k\) = cos a|Ai) + sina|A2), \k%) = 
— sin a|Ai) + cos a|A2). This is a two-dimensional subspace 
of the full Hilbert space supporting H. Here it has been cho- 
sen that {| Aii 2) defines the x axis on the Bloch sphere of 
FIGIi] hence Im(/c|A)=0. By confining | fci.2) to the span of 
|Ai 5 2/ then (using the fact it is diagonal in the measurement 
basis) one restricts interest to the component of <tg within this 
qubit space too: 



$e = e^-tc^ilfciX/cil + c^ fc 2|fc 2 )(fc 2 | + • • -> e 

/ c 2 e -i(\2-Xi)9 s 



-iH6 



s 2 -e~ l{ - x ' 2 ~ Xl ^sc 
C ^> k2 \ - e + l ( A 2-Ai)e) sc c 2 



where c (s) is cos a (sin a). We ignore elements of $0 
that project onto the remaining Hilbert space, orthogonal to 
{|Ai),|A 2 )}. Defining j3 = \ - (A 2 - \\)9, then angles 
{a, [3, 7} give an expectation value (if)o\H$$H\ipQ): 
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-4(X 1 -X 2 ) 2 s 2 [2a] s 2 [2 1 ]s 2 [fl 

(c[2(a - 7)] + c[2(a + 7)] + 2c \0\ s[2a]s[2 7 ] - 2)(c[2(a - 7)] + c[2(a + 7)] + 2c \0\ s[2a]s[2 7 ] + 2) 



(10) 



writing sin as V and cos as 'c'. J is optimized by angles 
{a, 7} 1 — ► 7r/4, independent of the value of f3 and giving a 
saturable bound: (ijj \H&gH\ipo) < (Ai - A 2 ) 2 . This is the 
upper bound on the Fisher information for any superposition 
of two eigenstates of the Hamiltonian. It is saturated by a 
probe state |V> t ) = (|Ai) + e l *\ X 2 ))/y/2, where X € [0, 2tt), 
see FIG[T] The result a 7r/4 dictates an optimal measure- 
ment scheme with components: 



(|Ai)±|A s 



(11) 



also in FIG[T] Other basis elements 6 {| k)} span an orthogo- 
nal subspace. It is significant that the optimal measurement is 
independent of (3, and hence the estimated parameter 9. 

Generalization to Higher Dimensions: The above result 
shows that for a given H, the maximal Fisher information is 
bounded from below by (Af — A^) 2 = ||-ff|| 2 where Af (Aj) 
is the max (min) eigenvalue of H, and \\H\\ = (Af — AjJ is 
the operator seminorm of the Hamiltonian J7). The variance 
has the seminorm as an upper bound: ||-ff|| 2 > 4A 2 H, cre- 
ating a bridge between the qubit result with that for a fixed 
\ifjo) in a higher dimensional space. We saw for a fixed |^o) 
that J < 4A 2 H. Therefore the qubit maximum variance 
state must be the universally optimal state over the full Hilbert 
space; it saturates the variance bound. Concisely: 



max 



h/>o),{|fc)} 



J 



\H\ 



(12) 



(This bound has been discussed previously in terms of quan- 
tum Fisher information Q.) A corollary of Eq.([T2]) is that no 
greater number of superposed energy eigenstates can be used 
as an input to improve on the Fisher information provided by 
the (qubit) maximum variance state ( | At ) + e lx | A^ ) ) / \/2. The 
optimal measurement set can be chosen as the one with two 
elements straddling the qubit subspace of extremal eigenval- 
ues in Eq.([TTj). Importantly, it retains independence from the 
true value of 9. 

Optimal Measurement Criterion: There may certainly exist 
more than one optimal measurement set for a given input, and 
we can define an optimal measurement as one that satisfies 
JC = (H) 2 , i.e. for measurements saturating J < 4A 2 H: 



A 2 ^ = K - (H) 



E 



r 2 Sk~ (v)c) 2 ^0 (13) 



Using Eq.(|9]l, for optimality <fik = (<p) c = — (if), Vfc. This is 
equivalent to a condition presented in [4|, lm(ipg\k) (k\ip±) — 
0, Vfc, where \4>±) = \ip) — (V'el'^IV'e) i s me component or- 
thogonal to \ipg) in the qubit space spanned by {\ipg), |^)}- 
Measurements that are superpositions of {\ipe), I^J-)}' e.g. 



those of [8], are generally only instantaneously optimal - as 
the state evolves, precision decreases. Demanding parameter 
independent optimality means <f>k = — (H) at all times, i.e. 



= -(H)e + a 



Vfc, 9 



(14) 



Enforcing this condition will limit the viable probes, measure- 
ments and dynamics. Starting with e l * fc (r + ircpk) = (k\ip) = 
—i(k\H\^) from just after Eq.(|2]i and substituting Eq.(fT4]i re- 
stricts the re-zeroed Hamiltonian, H = H — (H)I : 



(15) 



Here we have chosen an optimal |fc) basis, in which H has 
a real (imaginary) part denoted by R (I). Only the imagi- 
nary part determines dynamics |Q3), and \ipg) is confined to 
the null space of P. That H^ exp{-iH^ 66}\i) ) = 
implies [H^ R \ H^ 1 '] = 0. Therefore, in the optimal basis the 
Hamiltonian is block diagonal: H = H^ ® H^\ Note the 
null space of H^ needs to be at least two-dimensional for 
\ipo) to evolve at all - we have seen that maximum precision 
is possible in just such a qubit space. (Then if ( 7 ) must be pro- 
portional to the Pauli a y .) Complete measurements like M are 
not covariant; see ||5][T0l[T4| for a discussion of over-complete 
covariant measurements. 

Stability at Optimal Point: Fulfilling the conditions above 
may indeed give maximum precision but what if small devi- 
ations in the measurement orientation lead to a dramatic re- 
duction in the Fisher information? If the measurement basis 
is rotated slightly by |fc) 1— > exp{— ihduj}\k) then the over- 
lap becomes (fc|"0e) l_ ► (k\e lh5ul e~ lH89 \ipg), combining dy- 
namics for measurement drift and state evolution. The in- 
formation complement is now a function of two variables, 
JC(9,cu). At a turning point we have dgJC — d u JC = 0, 
but this doesn't indicate much. However, at the global min- 
imum (gm) for parameter-independent evolution Eq.([T4} ap- 
plies and therefore (p^ = —dg(H) = 0,Vfc. Here derivatives 
are dgJC — dgK, — 0, meaning that to enforce the optimal- 
ity condition for all 9 values confirms a zero curvature of J 
in the direction of evolution. Denoting derivatives of 9 by a 
dot and of ui by a dash, the other second order derivatives for 
parameter-independent evolution are: 



d 2 JC 



dudO 

d 2 JC 
duo 2 



=-2(h)y: 

gm k 

= m 1 ? 

gm 



Pk 



b 'k , 



-2(H) 



k 



p' k 4>'k 



(16) 



Probes returning (H) = produce a minimum in JC for de- 
viations Sui from the optimal measurements - then d^JC = 
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FIG. 2: For a spin Hamiltonian H i— > J y the maximum variance 
state |?/^o } = (| j, +j)y + exp{ix}|j, —j) y )/V2 (sometimes called a 
NOON state QiDlII]) yields greatest precision: J = || j y \\ 2 = 4j 2 . 
Eigen-equations are J 2 \j, m); = j(j + l)|j, m)i and J;|j,m) 4 = 
m\j, m)i for i € {x,y,z}. The optimal measurement set is not 
unique; one is given by Eq.^TTJ (with Ai,a *— ► ±j) but others include 
J x . Spherical surfaces (a) and (b) are coloured by precision: Fisher 
information for rotational drift |fe) i— > exp{ij z ui z } exp{i J y u) y }\k) 
of optimal measurements away from J x in (a) or from the set of 
Eq.ljTTJl in (6). Blue regions correspond to precision below the clas- 
sical bound (associated with J7max for n = 2j individual spin 1/2 
particles, J < 2j). Red regions indicate supra-classical precision 
2j < J < 4j 2 where quantum correlations contribute to the pre- 
cision 1 18 1. The clamshell structure of (a) is characteristic of the 
measurement landscape at all j with the red zones or "hotspots" 
numbering 4j around the equator, each decreasing in size as j in- 
creases. In (a) precision is highest and curvature is zero along the 
equator, showing that measurements cos^j^ + sin£J z are optimal 
and parameter independent for dynamics J y . However, curvature is 
non-zero along lines of longitude and dependent on drift angles ui ytZ 
so \ipo) may offer greatest precision but low stability. Precision con- 
tours in (c) mark the classical precision boundary for NOON states 
of spin j = 3/2 (dashed) and 21/2 (unbroken). The angle 5u± in- 
dicates how much transverse drift can be tolerated while maintaining 
supra-classical precision. 

and d^JC > 0. (Precision is optimal in both evolution and 
drift variables.) Ideally (H) = d^ g (f>k — 0, Vfc, then optimal 
measurements are perfectly stable, i.e. JC has a zero Hessian. 
Well, (H) = is easily obtained by a recalibration of the zero 
energy point, H i— > H without changing the physics of the 
system. However, the mixed second derivative of (f>k cannot 
be zero if drift dynamics rotate the measurement basis closer 
to the eigenbasis of H. (Parameter information becomes zero 
when measurements are H eigenstates.) The ideal scenario of 
Eqs.(fT5} produces Ti[A1H^] = and maintaining this under 
drift e~ Koh Me l implies an orthogonality condition: 

Tr(# (/) [M,h]\ = (17) 

Summary and Outlook: The formalism we have developed 
incorporates all aspects of quantum parameter estimation ex- 
plicitly; probe IV'o), dynamics H, and measurement, {|fc)}, 
clarifying how precision is determined by the interplay of all 
three. The information complement JC was introduced, a mea- 
surement dependent functional that undermines the Fisher in- 
formation. At the global optimum for parameter-independent 



measurement, phases = arg(fc|7/>#) vary linearly with the 
parameter 8 in proportion to their average energy, restricting 
dynamics to imaginary Hamiltonians in the optimal measure- 
ment basis. 

Greatest possible precision within the Hilbert space 
spanned by the Hamiltonian exists in the qubit subspace of 
the maximal variance input state and ultimate precision over 
all probes and measurements is completely defined by the ex- 
tremal energy eigenvalues. No additional dynamical structure 
is relevant, nor is the dimension of the Hilbert space. The most 
precise measurement is also parameter free. 

Stability of optimal measurements was quantified in terms 
of the curvature of JC in the vicinity of its global minimum. 
For dynamics h causing unitary drift of measurement ori- 
entation the curvature indicates whether the apparatus is of 
pragmatic utility, quantifying its immunity to alignment er- 
rors at the optimal setting. States with zero average energy 
(H) are associated with optimal parameter-independent mea- 
surements that become suboptimal given measurement drift, 
at a rate determined by the magnitude of mixed second order 
derivatives of phases <fik- See FIG |2] for an illustration of the 
precision terrain for spin states proposed previously for pa- 
rameter estimation. 

In future, it may prove fruitful to develop the approach pre- 
sented here to incorporate evolution governed by completely 
positive maps and generalised measurements. This work was 
carried out under a contract with Mission Critical Technolo- 
gies at NASA Ames Research Center. The author thanks 
Hugo Cable, Gen Kimura and Vadim Smelyanskiy for useful 
discussions 
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